Two opposite hysteresis curves in semiconductors with mobile dopants 
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Recent experimental researches on semiconductors with mobile dopants (SMD) have reported 
unconventional hysteretic current-voltage (I-V) curves, which form dynamically in either one of the 
two opposite directions, the counter-figure-eight and figure-eight ways. However the fundamental 
theory for the formation of the two directions is still absent, and this poses a major barrier for 
researches oriented to applications. Here, we introduce a theoretical model to explain the origin 
of the two directions, and find that the two ways originate from the spatial inhomogeneity of the 
dopant distribution in SMD. The counter-figure-eight (figure-eight) way of the hysteretic curve is 
obtained when dopants are driven from far from (near) the Schottky interface to the opposite side 
in the SMD. This finding indicates that the directions of hysteretic curve in SMD can be controlled. 

PACS numbers: 73.40.Sx,68.47.Fg,73.30.+y 
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The successful application of semiconductor devices 
over a wide range of electronic devices is attributable 
to their unique electrical properties, which are sensi- 
tive to the configuration of their dopants. Generally, 
dopants are assumed to be immobile. What happens if 
the dopants are mobile? Recent research on semiconduc- 
tors with mobile dopants (SMD) such as oxygen vacancies 
(O v ) [H indicates that SMD demonstrates its conduc- 
tance change due to the alteration of the dopant distri- 
bution in SMD Q. As an example, consider a metal(M)- 
semiconductor(S) Schottky contact (n-type) as shown in 
Fig.QJa). When a positive (negative) external bias V ex t is 
applied to the metal, the dopants are pushed away from 
(attracted toward) the metal, as a result, the dopant dis- 
tribution is non-uniformly altered. Because the interfa- 
cial band structure of a semiconductor is sensitive to the 
concentration of the dopants, different dopant distribu- 
tion causes different conductance state of the semicon- 
ductor. This conductance change is applicable to var- 
ious fields of science and engineering such as resistive 
switching phenoemena memristive system 3l,neuro- 
science [5(, and non- volatile memory devices 6H9|, and 
thus it has received great interests. 

Such conductance change results in a hysteretic 
current-voltage (I-V) curve. One conventional explana- 
tion for the hysteretic curve is as follows Q. When a 
positive (negative) bias V+ (V_) is applied, the donor 
concentration rid becomes low (high) near the Schottky 
interface. Then the Schottky barrier width w s b increases 
(descreases) because w s b oc l/y/n& thus, the con- 
ductance decreases (increases) as denoted by (T) (@) in 
Fig. [Tf b) . This directional hysteretic curve is called the 
counter-figure-eight (cF8) curve (ill ] and has been found 

in many materials [IS El- 

However, the opposite way of curve (Fig. QJc)), called 
the figure-eight (F8) curve, has also been reported in 
many literatures [y, 0, [3]. Moreover, the coexis- 
tence of both directions in a single sample has been re- 



ported [H, [U [l5| . To understand the origin of the two 
directions, several experiments and heuristic arguments 
have been presented. For example, Yang ct al. sug- 
gested that the cF8 and F8 curves are derived from the 
top and bottom Schottky interfaces, respectively. How- 
ever, Shibuya et al. [ll| hypothesized that the cF8 curve 
arises from O v movement through conducting filaments 
inside the sample, whereas the F8 curve has a purely elec- 
tronic origin. Subsequently, the same authors suggested 
that the cF8 and F8 curves originate from the respective 
inhomogencous (or filamentary) and homogeneous distri- 
butions of O v parallel to the interface . 

Despite these experimental results and heuristic argu- 
ments, the origin of the two ways of hysteretic I-V curves 
has not been elucidated theoretically yet. In this Let- 
ter, we theoretically demonstrate that the two ways of 
I-V hysteretic curves intrinsically appear in the SMD, 
resulting from the spatial inhomogencous distribution of 
dopants. When most dopants are located far from (near) 
the Schottky interface, the cF8 (F8) curve appears. 

We first introduce a simple theoretical model in one 
dimension to illustrate the mechanism of the two ways 
in SMD depending on the initial dopant distribution. In 
this model, the semiconductor is in contact with the met- 
als located at x = and x = I to form the Schottky and 
Ohmic interfaces (Fig. IHd)) with the boundary condi- 
tions Ec(x = 0) = Eq and Ec(x = I) = 0, respectively. 
We consider two different cases of initial dopant density 
distributions: dopants arc located 1) far from (Fig. Q] 
(d)) and 2) near (Fig. [T](c)) the Schottky interface. For 
simplicity, we assume that the dopant density distribu- 
tion pd(x) is constant in space. Then, for doped region, 
Pd = Q I (L — £) in the region [x — £, L] for the far-from- 
Schottky case and pd = Q/l in [x = 0, £] for the near- 
Schottky case, where Q is the total amount of dopants 
in a semiconductor and assumed to be a conserved quan- 
tity. For undoped region, pd = 0. This simplification is 
very useful to capture the essential mechanism of the two 
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FIG. 1. (Color online) (a) Diagrams of a O v -based SMD. 
Dopants can be repelled or attracted by applying a positive 
or negative bias, respectively, (b) and (c) are the schemat- 
ics for counter-Figure-8 and Figure-8 I-V hysteresis curves, 
respectively, (d)-(g) one-dimensional SMD model, (d) and 
(e) show the dopant density distribution pd{x) for the far- 
from-Schottky and near-Schottky cases, respectively, (f) and 
(g) show the effects of donor movement on the Schottky bar- 
rier for the far-from-Schottky and near-Schottky cases, re- 
spectively. 



ways of the hyteresis curves. We assume that the elec- 
trons are fully depleted in the doped region for analytic 
calculation. Non-constant Pd(x) case and not- fully de- 
pleted cases will be treated numerically later. Under this 
simplified circumstance, the position-dependent conduc- 
tion band Eq(x) can then be calculated by solving the 
Poisson's equation [lij], V 2 Ec{x) = ep sc (x)/e, where e is 
the electronic charge, p sc (x) is the space charge density, 
and e is the permittivity of the semiconductor. Note that 
p sc {x) = qpd(x), where q is the dopant charge. Here, we 
deal with the case q > (n-type semiconductor). 

We hrst consider the far-from-Schottky case. The Pois- 
son's equations for Ec(x) in the regions x < £ and 
x > £ become d 2 Ec{x)/dx 2 = and d 2 Ec(x)/dx 2 = 
qeQ / 'e{L — £), respectively Using the boundary condi- 
tions, Ec{x = 0) = Eq and Ec(L) = 0, and continuity 
at x = I, we can easily obtain Ec{x) in the whole range. 
Particularly for x < £, we obtain that 

E c (x)=f lX +E , whereA =-^-^p^. (1) 

Here, /i is the slope of Ec in the undopped region. If £ is 
initially located at © in Fig. [TJf) and a positive bias V+ 
is applied, £ increases as the direction ©— >©— >©. Then 
fi increases or the slope in the undopped region becomes 
less steeper (Eq. ([1])) as shown in Fig.QJf), which makes 
the Schottky barrier width w s h thicker. Therefore, the 



conductance decreases, which corresponds to the conduc- 
tance change denoted by ® in Fig. [lib). If a negative 
bias V*_ is applied to this low conductance state, I will 
change reversely as ©—>©->©. Then the conductance 
increases as denoted by © in Fig.QJb). This result agrees 
with the conventional explanation for cF8 curve. 

For the near-Schottky case, the calculation for Ec{x) 
can be performed similarly. The Poisson's equations for 
x < £ and x > £ become d 2 Ec(x)/dx 2 = qeQ/e£ and 
d 2 Ec(x)/dx 2 = 0, respectively. For x > £, we obtain 
that 

E c (x)=f 2 (x-L), where f 2 = -^ + ^. (2) 

If initial £ is located at @ in Fig. [T^g) , V+ makes £ in- 
crease as the direction @— >©^©. Then, by the similar 
explanation as the far-from-Schottky case, w s b becomes 
thinner as shown in Fig. [ljg) and the conductance in- 
creases, which corresponds to the conductance change 
denoted by © in Fig. []Tc) . If U_ is applied to this high 
conductance state, reverse process occurs, which causes 
the conductance decreases as denoted by © in Fig.[TJc). 
Therefore, this result verifies that F8 curve intrinsically 
appears in SMD without the assumption of the electronic 
function or the two Schottky interfaces. 

Using numerical simulation, now, we will show that the 
above analysis is still valid without the assumptions of 
non-uniform pd{x) and not-fully depleted doped region. 
For quantitative calculation, we use the parameters for 
Pt-SrTiC>3 contact. Let us consider a three-dimensional 
lattice (lattice constant a = 0.39 nm) whose lengths in 
x-, y-, and z-directions are L x , L y , and L z , respectively, 
as shown in Fig. HJa). Here, L x = L y = L z = 27.3 nm 
for manageable calculation. Two different metals are in 
contact with the lattice at x = and x = L x forming 
Schottky and Ohmic contacts, respectively. Donors (O v ) 
were distributed on the lattice depending on pd(x). Then 
the position-dependent conduction band Ec(xi,yj, Zk) 
can be calculated numerically by solving the Poisson's 
equation, V 2 E c (xi,yj, z k ) = ep sc {xi,y 3l z k )/e. How- 
ever, the calculation of Eq is not straightforward because 
Psc{x) qpd{x). Therefore, we use the self-consistent re- 
laxation method to obtain p sc and Ec simultaneously; 
we divide p sc into two parts: p sc = p + — p_, where 
p + and p_ are the densities of positive and negative 
charges, respectively. When the density of donors is 
sufficiently high compared with the hole density, p+ ~ 
qpd(x). p-{xi,yj, z k ) corresponds to the density of elec- 
trons and is determined by the following equation (To| : 
p-(xi, yj ,z k ) = 27V C /0F J^dr) r?°- 5 /(l + exp [ v -p{E F - 
Ec(xi, yj, z k )}}), where /3 is the inverse temperature and 
N c is the effective density of the states in the conduc- 
tion band. We used N c = 2.5 x 10 19 cm -3 because 



10 cm for many semiconductors [If 
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also confirmed that the essential feature of the simula- 
tion is not changed by variation of N c . Note that we 
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FIG. 2. (Color online) (a) Configuration of the simulation. 
The three-dimensional lattice (semiconductor) is in contact 
with different metals to form a Schottky contact at x = 
and an Ohmic contact at x = L x . Donors are denoted by 
yellow circles, (b) Periodic potential energy of the donors 
(grey curve), where local minima correspond to the lattice 
sites. Uo is the energy barrier height against the movement 
of a donor. When a negative (positive) voltage is applied, the 
periodic potential energy increases (decreases), as indicated 
by the orange (purple) curve. The slope E corresponds to the 
applied electric field. Donors move according to the hopping 
probabilities po, p+i, and p_i (Eq. (|4|). 



set Ep = V ex t hi the semiconductor and assume that the 
barrier height at the M-S interface is independent of the 
dopant density [17j | ■ Thus, we can set up the bound- 
ary conditions at x = (ideal Schottky) and x = L x 
(ideal Ohmic) interfaces as Ec{0,y,z) = 0.9 eV [l8| and 
Ec{L x ,y, z) = T4 x t, respectively. Here, we neglect the 
image-charge effect for the sake of simplicity. Inserting 
p+{xi,yj, Zk) and p-(xi,yj, z k ) into the Poisson's equa- 
tion, we obtain E c (xi,yj, z k ) and p sc (x l ,y J , z k ) simul- 
taneously To confirm the validity of this technique, 
we calculate Ec{xi,yj, z k ) for a silicon semiconductor 
with various doping concentrations. The results are pre- 
sented in Supplemental Material (SM) 1. The resulting 
Ec{xi,yj, Zk) describes the conduction band which is dis- 
torted in the insulating region of the Schottky barrier. 

Next, using the obtained Ec(xi,yj,z k ), the electric 
current / of the major carriers (i.e., electrons) can be 
estimated with the following formula [lj| : 
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l + exp{/3(Z-E x -V ext )) 



where m e is the free electron mass, h is Planck's con- 
stant, and £ = max(_Ep — Ec)- Pj, k (E x ) is the transition 
probability that an electron with x-directional energy E x 
will tunnel through the Schottky barrier at y = yj and 
z = Zk- In the discrete lattice, Pj k (E x ) can be written as 



Pj,k{E x ) « eyLp^-aJ2 l a^E c {x i ,y j ,z k ) - E x j , where 
the summation index i extends over all cases satisfying 

E c (xi,yj,z k ) > E x and a = 1.025 cV~ - 5 A~\ 

We assume a simple hopping motion along the in- 
direction for the donors under a periodic potential with 
a barrier height Uo , as shown in Fig. [U[b) , because it is 
widely accepted that O v must overcome an energy bar- 
rier to move to the next lattice site [20[ . We also assume 
that a constant electric field E = —V cxt /L x is formed 
throughout the semiconductor when V cx t is applied be- 
tween two electrodes. The validity of the constant E- 
field approximation is discussed in SM2. Thus, when a 
negative (positive) V ey ± is applied, the periodic poten- 
tial energy for the donors increases (decreases) with a 
slope of E, as shown in Fig. [5{b). Then, the heights 
of the left and right energy barriers, compared to the 
local minimum, become approximately Uq — aE/2 and 
Uo + aE/2, respectively. The probability of remaining at 
the original site Xi (po) is given by the probability that 
the donor cannot overcome a lower barrier among the 
two. So, po = l-exp(-/3(£/ -a|.E|/2)). When V eyct > 0, 
the probability of moving to site (p_i) is the half 
of the probability that the donor overcomes the left or 
higher barrier (another half of the probability should be 
counted for moving to the opposite direction). So, p_i = 
0.5exp(-/3([/ + a|£;|/2)). Then p +1 = 1 - p - P-i = 
exp(-/3(U - a\E\/2)) - 0.5 cxp(-/3([/ + a\E\/2)). Simi- 
larly when Vcxt < 0, p +1 = 0.5 exp(-/3(C/ +a|i?|/2)) and 
p_i = ew(-P(U -a\E\/2))-0.5exp(-P(U + a\E\/2)). 
By combining these two cases, we obtain 

p+i = 0.5e- f3Uo \ e ^ E W 2 + 2sgn(V r xt)sinh(/3|£:|a/2; 



= 0.5e 



-PU 



,P\E\a/2 



2sgnG4xt)sinh(/3|£|a/2) 



(3) 



P-i 

po = 1 _ e -P(Uo-\E\a/2) t (4) 

where sgn(x) = —1, 0, and 1 when x < 0, x = 0, and 
x > 0, respectively. For simplicity, we consider only a 
hardcore repulsion interaction between the two donors. 
Here, we adopt the thermal acceleration mechanism ap- 
plied for SrTiC>3 2l|. So, high temperature ft ~ 15 cV _1 
(800 K) can be used for our simulation with Uq = 1.01 
eV [2l|. From Eq. (U]) it is obvious that transformations 
as /? — > P/b, Uo — > bUo, and E — > bE do not change the 
hopping probabilities. Thus, we use ft = 40 eV -1 (300 K) 
and U = 0.379 eV instead of P = 15 eV" 1 and U = 1.01 
eV. Here, the attempt frequency for the hopping is 10 13 
Hz Q. 

Using the above equations, the case in which the 
donors move from the Ohmic to the Schottky interface 
can be simulated. Initially, the donors were uniformly 
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FIG. 4. (Color online) If the donors are initially distributed 
in the far-from-Schottky region (a), the voltage sweep results 
in a cF8 directional G-V curve (b). If the donors are initially 
distributed in the near-Schottky region (c), a F8 direction 
(d) is obtained. By applying a large negative bias to the 
lattice exhibiting a F8 direction, we can attract donors into 
the near-Schottky region, and then the direction of G-V curve 
will change to F8 way. The opposite effect can be obtained 
by applying a large positive bias. 



FIG. 3. (Color online) (a) Changes in the donor density dis- 
tribution from t = to 14 /is when a negative bias is applied, 
(b) Changes in the conductance as a function of time, (c) and 
(d) show changes of the Schottky barrier when most donors 
are distributed in the far-from-Schottky and near-Schottky 
regions, respectively. In all figures, red, cyan, blue, green, 
and gold are used to represent data collected at t = 0, 1, 7.5, 
10, and 12 /is, respectively. 

distributed with a density of 10 19 /cm 3 [2l[. Using 
Eq. ((4]), we pushed the donors toward the Ohmic inter- 
face by applying a positive bias, the red curve in Fig.^a). 
Then, we applied a negative voltage V ext = —1.875 V to 
attract donors towards the Schottky interface. Here, the 
donor density at xt is defined as ~p{xi) = n(xi)/(L y L z ), 
where n(xi) is the number of donors at the x = Xi plane. 
Here, e = 100en in high electric field (en the permit- 
tivity in free space), with periodic boundary conditions 
in the y- and z-directions. Fig. [3Ja) shows the time- 
dependent distribution of the donors. The distribution 
moved toward the Schottky interface over time. 

The conductance G{= I/V cx t) during the attraction 
process is calculated at 0.1 V as a function of time t. 
As indicated in Fig. |3^b), the G-t plot can be divided 
into two regions: for t < 7.5 /is, G increases as a func- 
tion of t, and for t > 7.5 /is, G decreases. When com- 
paring the distributions shown in Fig. [3Ja), G increased 
(decreased) when most dopants were distributed in the 
far-from-Schottky (near-Schottky) region. 

The two different t-dependences of R come from the 
different modulation behaviors of the Schottky barrier 
during the attranction process. For each t, we obtained 
Ec{xi) by calculating Ec{xi,yj, Zk) at V ex t = and aver- 
aging over yj and Zk- Fig.JHJc) shows Ec{xi) when most 
dopants were distributed in the far-from-Schottky region 



(i.e., t < 7.5 /xs). In this case, the pulling of the donors 
toward the Schottky interface resulted in a decrease in 
the Schottky barrier width, and G increased. Fig. (3][d) 
presents the case Ec{xi) where most dopants were dis- 
tributed in the near-Schottky region (i.e., t > 7.5 /is). In 
this case, the attraction of the donor increased the Schot- 
tky barrier width rather than decreasing it. These results 
agree with those of the one-dimensional SMD model. 

We also simulated G-V curves under a repetitive volt- 
age sweep, with different initial donor distributions. 
Here, it took 0.1 fis for each voltage point and voltage 
gap is 0.027 V. During the voltage sweep, G is calculated 
at 0.1 V. When most donors were initially distributed in 
the far-from-Schottky (near-Schottky) region as shown in 
Fig. 0Ja) (Fig. Hfc)), G decreases (increases) at positive 
bias and increases (decreases) at negative bias as shown 
in Fig- HJb) (Fig. Hfd)), which corresponds to the direc- 
tion of Fig. [Ub) (Fig. [He)) ( see a ^ so movie clips). Fur- 
thermore, the direction of the G-V curve can be changed 
by applying a large bias, as indicated in Fig. |4l 

In conclusion, we introduced the SMD model which 
demonstrates that two opposite hysteresis curves intrin- 
sically appear in the SMD due to the inhomogeneous 
dopant density distriubtion. From this theoretical anal- 
ysis, we can control the type of the I-V curve by mod- 
ulating the mobile dopant distribution. The theoretical 
result we obtained in the Letter may become a funda- 
mental basis for further development of SMD. 
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